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The Crystal-GRID method is used to study 
interatomic collisions at low energy in 
metals and such to probe the repulsive in- 
teratomic potential. Line shapes of 
gamma rays, emitted by the recoiling 5, Ni 
isotope after thermal neutron capture in 
Ni single crystals, were measured and com- 
pared to results obtained by molecular 
dynamics simulations of the slowing down. 
The same procedure is also used for re- 
coiling 57 Fe and 54 Cr atoms in Fe and Cr 
single crystals, respectively. Different po- 
tentials (including several from the embed- 
ded atom method) are investigated using 
the observed fine structure of the line shape 
which depends on the crystal orientations. 



From the detailed study of the lineshapes 
measured in two different orientations, a 
new potential is then derived for each ele- 
ment. Nuclear state lifetimes for the ex- 
cited isotopes are also deduced with a 
higher precision than obtained with stan- 
dard nuclear techniques. 
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1. Introduction 



The gamma-ray induced doppler broadening (GRID) 
method [1] is based on the observation of Doppler shifts 
produced by the motion of radioactive nuclei in solid 
state targets. The motion of atoms is produced by de-ex- 
citation following thermal neutron capture. Usually after 
this type of reaction, the newly formed isotope will 
decay by emission of gamma-ray cascades down to the 
ground state. The GRID technique involves the measure- 
ment of the energy of the second gamma ray in a two 
gamma cascade. The first gamma ray induces a recoil to 
the nucleus in the opposite direction of emission. If the 
nucleus emits a second gamma ray while still in flight, 
the measured energy of the latter will be Doppler 
shifted. As the gamma rays are emitted isotropically, this 
results in a Doppler broadened line shape rather than a 
line that is Doppler shifted. This small Doppler broaden- 



ing can be measured by the two-axis flat crystals spec- 
trometer GAMS4 [2] installed at the Institut Laue 
Langevin (ILL) in Grenoble (France) in a ILL/NIST 
collaboration. Typical values for the recoiling atom en- 
ergy analyzed with the transition metals are in the order 
of 200 eV to 400 eV. With this kinetic energy, the atom 
has enough energy to move from its lattice position and 
make collisions with its neighboring atoms. Before the 
excited nuclei emits the second gamma ray (the proba- 
bility of emission is determined by the nuclear state 
lifetime), the atom will slow down due to interactions 
with its neighbors. The slowing down process is simu- 
lated by molecular dynamics programs where different 
interatomic potentials can be implemented. The result- 
ing Doppler line shape is compared to the measured one 
via a fitting routine [3]. The Crystal-GRID technique [4] 
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uses single crystals as targets, and was developed to 
study and to construct new interatomic potentials. Here 
we apply this method to analyze three transition metals: 
nickel, iron, and chromium. 

The Doppler broadened line shape depends on three 
contributions: i) the lifetime of the nuclear state, which 
gives the probability of emission of the second gamma 
ray at a given time, ii) the orientation of the crystal with 
respect to the direction of observation, which will be 
discussed in Sec. 2 and iii) the slowing down process, 
which will be discussed in detail in Sec. 3. The results 
of the different slowing down theories for the three 
transitions metals are summarized in Sec. 4. Section 5 
deals with the construction of a new potential derived 
from the measured Crystal-GRID data. 



2. Crystal-GRID Technique 

The gamma ray induced doppler broadening mea- 
surements are realized with single crystals as target. The 
targets consists of three crystals with a dimension of 
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2 mm X 18 mm X 20 mm having a total mass of 10 g 
to 15 g. The targets are placed oriented inside the nu- 
clear reactor and are irradiated with a neutron flux of 
5 X 10 14 cm" 2 s" 1 . The different isotopes and nuclear 
reaction studied in this experiment are listed below: 
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The process involved in the production of a Doppler 
broadened gamma-ray line shape are described in Fig. 1 . 
First (i), neutron capture forms an excited nucleus. The 
newly formed atom will deexcite by emitting a first 
gamma ray (ii). This emission will induce a recoil to the 
nucleus and the atom starts to move in the crystal (iii). 
If a second gamma ray is emitted during the flight (iv) 
at a time depending on the nuclear state lifetime (t), it 
will be Doppler shifted due to the emitting atom veloc- 
ity. The GAMS4 spectrometer measures the energy of 
the second gamma ray with respect to a fixed direction 
and records this Doppler shift leading to the line shape. 
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Fig. 1. Description of the nuclear reaction and the nuclear state levels and transi- 
tions of interest for the Ni isotope (right) and the associated atomic events (left). 
For simplicity, the motion of the recoiling atom in the crystal is restricted in one 
dimension. 
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If the recoil is induced by only one primary gamma ray 
(no cascade or side feeding) then the kinetic energy (E r ) 
of the atom at the beginning of the recoil is given by: 



E r = E^llMc 2 



(2) 



where E yl is the energy of the first gamma ray emitted 
and M is the mass of the recoiling atom. Table 1 shows 
the nuclear levels as well as the different transitions 
characteristics investigated during this analysis. All of 
the nuclear levels used are primary almost 100 % di- 
rectly fed, which gives rise to a unique recoil velocity. 
As the nuclear state lifetime can be known by other 
means (DSAM, etc.), and as the direction of observation 
with respect to the orientation of the crystal is fixed, the 
data analysis gives valuable information on the slowing 
down process. 

Specially adapted molecular dynamics (MD) pro- 
grams [5] calculate the trajectories and velocities of the 
recoiling atom depending on the interatomic potential 
by solving the equations of motion of a set of atoms. 
Once the trajectories of the recoiling atoms have been 
simulated, the only free parameter left over, which can 
be varied during the fitting procedure, is the lifetime of 
the nuclear state. If the lifetime is very short, the atom 
will have a great velocity when it emits the second 
gamma ray as no collision has taken placed and the atom 
did not have sufficient time to slow down. Therefore the 
Doppler line shape will be more broadened. 

To reconstruct the line shape from the trajectories and 
velocities given by the MD simulation, the scalar 
product between the velocity and the direction of obser- 
vation (given by the direction of the spectrometer) is 
calculated. This gives the Doppler shifted energy of the 
emitted gamma ray. Then a summation over all trajecto- 
ries is realized weighted by the number of gamma rays 
emitted at that particular time. The probability of 
gamma emission is given by the exponential decay law 
assuming a lifetime value. Figure 2 shows the simulated 



line shape obtained for different lifetime values and 
different crystal orientations. 

Once the potential has been chosen, the trajectories 
and velocities of the recoiling atoms are stored, the only 
free parameter is the lifetime. This value is optimized in 
order to reproduce the measured Doppler broadened 
line shape. This fitting procedure returns a value for the 
lifetime and a x 2 P er degree of freedom which tells 
about the agreement/disagreement between the mea- 
sured and simulated line shape. A perfect match be- 
tween the two line shapes will result in a \ 2 equal to 
one. 

The crystal orientation influences the Doppler broad- 
ened line structure as blocking and channeling of the 
recoiling atom due to the ordered structure of the atoms 
in the target influences the trajectories and the velocities 
of the recoiling atoms. Therefore the line shape will 
show different patterns for different orientations. Ni, Fe, 
and Cr have simple fee or bec crystalline structure. Two 
orientations for each transitions metals were investi- 
gated, and are illustrated in Fig. 3. These two orienta- 
tions were chosen based on the crystalline structure, 
because they should show the biggest variation in the 
Doppler broadened line shape, as the distance separat- 
ing two closest neighbors is quite different for these 
orientations. 



3. Interatomic Potentials 

As the line structure and the fitted lifetime depend 
strongly on how the atom moves inside the crystal, the 
interatomic potential plays a great role in the data anal- 
ysis. There exists in the literature many interatomic po- 
tentials for metals and many ways to elaborate them. 
They are usually derived from bulk properties and fitted 
to experimental data when measurable. Many parame- 
trised forms and formulas (exponential, polynomes, 
spline knots, etc.) were developed in order to reproduce 



Table 1. The reaction and levels used for the three transitions metals 



Reaction 



(barn) 



A 

(%) 



(keV) 



E y > 
(keV) 



(%) 



X/ in /2/ 0l 

(%) 



(is) 



58 Ni(n,7) 59 Ni 4.6 68.3 4140.34 

2414.97 



56 Fe(n,-/) 57 Fe 2.63 91.7 3427.67 

1725.38 



3719.75 



3675.23 


0.96 


100 


6.5 (1.4) 


1950.05 


1.7 


96 


53.4 (7.2) 


2721.17 


1.37 


100 


4.33 (0.87) 


1725.29 


6.3 


86 


47.6 (7.2) 


2239.07 


13.9 


98 


<25 


3719.75 


5.04 


97 


<43 



' Reference [6] for Ni, [7] for Fe, and [8] for Cr. 
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Fig. 2. Simulated line shape obtained with the ZBL potential for different lifetime values and different crystal orientations for Fe. 
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Fig. 3. Direction of observation for the two crystalline structures. Ni 
has a fee structure shown on the left with the two orientations mea- 
sured and Fe and Cr have a bee crystalline structure. 

the value of the pair potential as a function of inter- 
atomic distance. 

Table 2 gives a list of some of the potentials used for 
the slowing down of atoms in metals. Only the best four 
potentials investigated are given in the table. These po- 
tentials are derived so to reproduce different characteris- 
tics of bulk properties, such as cohesive energy, lattice 

Table 2. Interatomic potential references 



Potential name 

Ni interaction 

BM 

ZBL 

EAMVC 

EAMFBD 

Fe interaction 

BM 

ZBL 

EAMVC 

EAMGA 

Cr interaction 

BM 

ZBL 

EAMGA 

EAMWB 



Authors 



Born Mayer 

Ziegler, Biersack, Littmark 

Voter and Chen 

Foiles, Baskes, Daw 



Born Mayer 

Ziegler, Biersack, Littmark 

Voter and Chen 

Guellil and Adams 



Born Mayer 

Ziegler, Biersack, Littmark 
Guellil and Adams 
Wang and Boerker 



Ref. 



[12] 
[13] 

[141 
[151 



[8] 
[9] 
[10] 
[16] 



[8] 
[9] 
[121 
[171 



constant, elastic constants, properties of stable struc- 
ture, etc. All of the mentioned potentials are pair poten- 
tials and are only dependent on the interatomic distance, 
even the recent developed embedded atom method 
(EAM) potentials [9]. In this new approach developed 
for metals, the atom is seen as embedded in the electron 
density caused by the neighboring atom and the total 
energy of the system can be described by the summa- 
tion over all atoms of an embedding function evaluated 
for an electron density and a core repulsive term. An 
approximation to have the force acting on one atom is 
needed in order to avoid the double summation (sum- 
mation over all electron density and summation over all 
atoms). The approximation replaces the electron density 
by an average electron density. A Taylor expansion 
over this average electron density is then realized for 
the evaluation of the needed parameters [10]. Figure 4 
shows the form of the different potentials for the Ni-Ni 
interaction. More potentials are given in the figure in 
order to illustrate the wide range of potentials found in 
the literature. A detailed presentation of all the poten- 
tials can be found in Ref. [11]. 



4. Results 

Many of the mentioned potentials can reproduce 
quite well bulk properties as they were designed for this 
purpose, especially the EAM potentials, however not all 
of them are suited for the slowing down in the energy 
range of a few hundreds of eV. The Crystal-GRID tech- 
nique can make a selection of the best interatomic po- 
tential for metals by comparing the lifetime value ob- 
tained and the \ 2 P er degree of freedom which stands 
for the agreement/disagreement between the stimulated 
and measured Doppler broadened line shape. The result 
for the different slowing down theories are listed in 
Table 3 for the three transitions metals. 



Table 3. Lifetime results for different potentials obtained with Ni, Fe, and Cr crystals 





E, 


Known lifetime 


BM 


ZBL 


EAMFBD 


EAMVC 




(keV) 


(fs) 


(fs) 


(fs) 


(fs) 


(fs) 


Nickel 5 'Ni 


4140.3 


6.5(1.4) 


4.21(0.39) 


6.88(0.47) 


11.89(0.96) 


7.95(0.55) 




2414.9 


53.4(7.2) 


26.24(0.95) 


46.30(1.60) 


94.18(4.01) 


53.06(1.80) 




x 2a 




1.1384 


1.1371 


1.1374 


1.1371 


Iron 57 Fe 


3427.7 b 


4.33(0.87) 


1.95(0.14) 


2.82(0.22) 


3.58(0.30) 


3.96(0.34) 




1725.4 


47.6(7.2) 


18.30(0.52) 


30.92(0.84) 


38.34(1.08) 


52.13(1.39) 




x 2a 




1.1298 


1.1292 


1.1301 


1.1288 


Chromium 54 Cr 


3074.0 


<24.5 


5.52(0.15) 


8.30(0.24) 


7.57(0.24) 


14.72(0.53) 




3719.9 


<43.3 


16.76(1.49) 


26.37(2.43) 


27.09(2.62) 


44.48(4.77) 




x 2a 




1.0854 


1.0859 


1.0918 


1.0817 



* X" taken for the shortest nuclear state level. 

b Gamma-gamma correlation was neglected in this analysis, despite the fact that it might occurs for this very short lived level, due 

to a lack of experimental knowledge of the mixing ratios. 
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Fig. 4. Theoretical interatomic potentials for the Ni-Ni interaction. All the potentials are placed in the figure from the most (BM) to the least 
(EAMATVF) repulsive potential. For the name of the potentials refer to Table 2. 



For nickel and iron, the best potential, which repro- 
duces the known nuclear lifetime and which gives the 
smallest \ 2 P er degree of freedom, is the embedded 
atom method potential from Voter and Chen (EAMVC). 
For chromium, the best candidate is the embedded atom 
method potential developed by Wang and Boerker 
(EAMWB). Figure 5 shows the lifetime value obtained 
with the different potentials and the known lifetime 
found in the literature for the three transition metals. 
Among the list of interatomic potentials found, some 
give lifetime value 3-4 times higher or lower than the 
known lifetime. This shows that these potentials are not 
suited in this energy domain. 

The orientation of the crystal also influences the 
Doppler broadened line structure. Figure 6 clearly 
demonstrates that the blocking and channeling due to 
the ordered atom position in a crystal gives rise to differ- 
ent line profiles. For example, in a bcc crystal (in the 
case of chromium), the atom is more free to move in the 
[110] direction, because the closest neighbor is further 



(see Fig. 3) than in the [100] direction. Therefore, the 
line profile will be more broadened in this particular 
direction compared to the [100] orientation profile. 



5. Construction of a New Potential 

From the GAMS4 data a new potential can be 
derived. This potential is based on the known potential 
given by Ziegler et al. This ZBL potential is a Coulomb 
screened pair potential with several fixed parameters. 
The potential is given by 



V(r»): 



Z,Z, 



<P(r,j) 



4tt£ r,j 
where the screening function <P has the form 

<P(r tj ) = X c, exp(^ 
1=1 \ a s 



(3) 



(4) 
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Fig. 5. Fitted lifetimes for the three transition metals. Two nuclear state levels were analyzed. The dotted line represents the known lifetime range 
found in the literature. 
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and the screening length a s is defined by 
a s = 0.8853a B (Z1 + Zff 



(5) 



with a B = 0.529 A the Bohr radius and Z tJ the atomic 
number of the atoms i andy respectively. All the parame- 
ters ici, di, x, y) for the evaluation of the ZBL potential 
are listed in Ref. [12]. As we are dealing with atoms 
recoiling in the same material, the x and y parameters 
can be considered to be one unique parameter. For the 
construction of the new potential, all the parameters are 
kept identical, except the x parameter, which is varied 
from 0.05 to 0.40. The original x value is equal to 0.23 
[9]. An optimization of the x parameter is realized by 
finding the minimum for the \ 2 P er degree of freedom 
obtained by comparing the simulated with the measured 
line shape. Figure 7 shows the behavior of the \ 2 as a 
function of the x parameter for the three transition 
metals. The best x value for the nickel and iron crystal 
is equal to 0.26 and for chromium to 0.31. The new 
derived potential has a form similar to the best slowing 
down theories found from the analysis when the selec- 
tion was realized with the lifetime and \ 2 P er degree of 
freedom as criteria. The construction of the new poten- 
tial from the Crystal-GRID measurement is independent 
of the nuclear state lifetime which in some case is not 
known with sufficient precision by other methods. 



6. Conclusion 

We have investigated the interatomic potentials for 
three transition metals using the Crystal-GRID tech- 
nique. This new method allows one to make a selection 
of the best slowing down description among a set of 
interatomic potentials and also permits to construct new 
potentials. Two criteria are preponderant for this analy- 
sis: i) the nuclear state lifetime and the \ 2 P er degree of 
freedom which is an indication of the agreement or 
disagreement between the simulated and measured line 
shape. Thus of a multitude of interatomic potentials for 
metals, the Crystal-GRID method shows that not all of 
them are valuable in the energy domain up to a few 
hundreds of eV. For this particular energy domain, 
where very few experiments exist, the Crystal-GRID 
method gives an opportunity to test the veracity of the 
slowing down of atoms in metals. This nuclear tech- 
nique based on the observation of Doppler broadening 
resulting from gamma emission after neutron capture 
also gives lifetime value for 6 nuclear levels with a much 
higher precision than standard nuclear techniques. 
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